Hydrodynamic description of neutrino gas 



Levan N. Tsintsadze 
Venture Business Laboratory, Hiroshima University, Higashi-Hiroshima, Japan 

Abstract 

The system of neutrino-antineutrino (z/P) - plasma is considered taking into 
account their weak Fermi interaction. New fluid instabilities driven by strong 
neutrino flux in a plasma are found. The nonlinear stationary as well as 
nonstationary waves in the neutrino gas are discussed. It is shown that a 
bunch of neutrinos, drifting with a constant velocity across a homogeneous 
plasma, can also induce emission of lower energy neutrinos due to scattering, 
i.e. the decay of a heavy neutrino vh into a heavy and a light neutrino i^l 
[vh — > vh^l) in a plasma. Furthermore we find that the neutrino production 
in stars does not lead in general to energy losses from the neutron stars. 



I. INTRODUCTION 

Our understanding of the properties of neutrinos in a plasma has recently undergone 
some appreciable theoretical progress [1-5]. The interaction of neutrinos with a plasma 
particles, the creation of vv pairs, the emission of neutrinos due to the collapse of a star are 
of primary interests [1,6-9] in the description of some astrophysical events such as supernova 
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explosion. One of the key processes upon the explosions are the large-scale hydrodynamic 
instabilities as well as v driven plasma instabilities. These processes are also believed to have 
occurred during the lepton stage of the early universe. During the formation of a neutron 
star the collapsed core of the supernova is so dense and hot that v and V are trapped and 
are thus unable to leave the core region of the neutron star. The rates of escape of v and 
p are very small, inside the star an equilibrium state is reached, which includes the 
concentration. Recently, the neutrino transport phenomenon during the Kelvin-Helmholtz 
phase of birth of a neutron star in the diffusion approximation was investigated by Pons 
et al. [10]. A detailed knowledge of transport properties of the neutrinos in these extreme 
environments must include the neutrino-electron Fermi weak interaction coupling. There 
are different mechanisms involved in the creation of the i^, as well as pairs. The collective 
effects of the stellar plasmas can significantly alter the production rate of neutrinos. It is 
widely thought that the most dramatic plasma process is the decay of photons as plasmons 
into neutrino pairs. It was pointed out by Adams et al. [11] that the neutrino pair radiation 
can be the dominant energy loss mechanism for the plasma of the very dense stars, as well as 
white dwarfs, red giants and supernovas. It has been also shown [12] that a vh can undergo 
radiative decay into a photon and a {i^h I^l) in the presence of a strong magnetic field, 
with the strength greater than the critical value given by Her = fn^c^/eh = 4.41 • 10^^ Gauss. 
The increase (decrease) of the u (p) energy within the plasma was also demonstrated [4]. 
This may have a significant and potentially detectable effect. Observationally, the recent 
results from the Kamiokande group provide strong evidence for the existence of neutrino 
oscillations. These results, together with the increased confidence in abihty to produce 
and manipulate intense muon beams makes feasible the future neutrino factories based on 
muon storage rings. A muon storage ring, as presently envisioned, would have energy in 
the 10-50 GeV range, and produce directed beams of intense neutrinos, and naturally such 
an intense stream of neutrinos can exist in astrophysical and cosmological plasmas. Thus, 
it is of interest to examine if new physical processes that are caused by the intense flux 
of neutrinos, can appear in such a plasma, i.e., of particular conceptual interest are effects 

2 



which have no counterpart in vacuum. 

In this paper we consider a medium of neutrino gas and plasma in semitransparent 
regions, where collective process plays an important role. In our consideration we assume 
that the interaction between neutrinos is weak in comparison with the neutrino-electron 
interaction, i.e., the neutrino gas is ideal. We can also ignore the spin of neutrinos in 
an isotropic plasma, since in this case the energy of neutrinos is independent of the spin 
operator. With this modeling, here we show that there exists the opposite effect that has 
been proposed by others [2,11]: namely that the neutrinos production does not generally 
lead to energy losses from a hot and dense system. 



II. BASIC EQUATIONS 

In the following we shall demonstrate the excitation of longitudinal oscillations in neu- 
trino flux-plasma medium. For this purpose, we make use of hydrodynamic equations for 
the neutrino gas, which have been derived by Tsintsadze et al. [4]. These equations for the 
neutrino density n,^ and the fluid velocity Ui, take the following forms 



dt 



+ div(n^u,y) — (1) 



at h < Uk > noe 

where U = \/2Gpnoe describes the universal Fermi weak interaction, Gp = I0~^^erg.cm^ is 
the Fermi coupling constant, — I (i^^' electron plasma density. 

The electron continuity equation takes the form: 

+ div{neVe) = 0. (3) 

For the equation of motion of electrons in the presence of a neutrino gas, we have 

^ = eV0 - V2GFVn, , (4) 

where we have neglected the pressure term in comparison with the first term on the right 
hand side under assumption that rjj^/L"^ <^ 1, where roe is the Debye length of the electrons. 



and L is the characteristic length of the system. Note that the last term in Eq.(4) is the 
neutrino ponderomotive force. A more elaborate discussion on the neutrino ponderomotive 
force is given in [4]. 

The electrostatic potential is determined from Poisson equation 

= 47re(ne - rij) . (5) 

We first consider a situation in which the ions are in an equilibrium state, i.e. rii — noi, 
where noi is the equilibrium value of the ion density. In this consideration the range of 
frequency of small perturbation of quantities of the medium is much less than the Langmuir 
frequency of the electrons. Thus, we can neglect the inertial term in Eq.(4), and define (f) 
from it. Finally, we obtain from the Poisson equation the expression for the electron density 

— = 1 + -; — 5 V , (6) 

noe 47re^noe uq^ 

where uqi, is the equilibrium value of the neutrino density. Substituting expression (6) into 
Eq. (2) we get the following equation, which except quantities of the neutrinos contains only 
the electron charge 

^ + (C-VR = -/3V(W) , (7) 
ot 

where = G%c'^nQy/27re'^h <ujk>- 

Equations (1) and (7) form a closed set of equations that describes the hnear as well as 
nonlinear waves in the neutrino gas. 



III. CHERENKOV TYPE OF EMISSION OF LIGHT NEUTRINOS 

Let us consider the propagation of small perturbations in a homogeneous plasma- neutrino 
beam. To this end, we linearize Eqs.(l) and (7) with respect to the perturbations, which 
are represented as — Uq^ + dn^, Un — Uq^ + 5ui,, where the suffix denotes the constant 
equilibrium value, and Srii,, 5ui, are small variations in the wave. After linearization of 
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Eqs.(l) and (7), we will seek plane wave solutions proportional to expi{qf — D,t). We can 
then derive the following dispersion relation 

(n - quo^f = -(3q' . (8) 

This leads to a solution of the form Q, — ReO, + ilmQ such that: 

ReQ, — quoi^cosQ (9) 

and the growth rate for the unstable branch becomes: 




This solution clearly describes the emission of the low- frequency neutrinos inside a resonance 
cone (cos^ = ) similar to the well-known Cherenkov emission of electromagnetic waves 
by charged particles moving in uniform medium with a velocity larger than the phase veloc- 
ities of the emitted waves. On the other hand, strictly speaking the effect considered here 
is physically quite different from the usual Cherenkov emission by charged particles. The 
point is that the neutrino flux can no longer excite Langmuir waves and energetic (heavy) 
neutrino scatters of the electron density perturbation (caused by the ponderomotive force of 
neutrinos) and emits low (light) energy neutrino. We specifically note here that such decay 
process vh ^h^l is forbidden without a plasma. 

Let us estimate the growth rate (10) for some typical neutrino parameters. Accordingly, 
wc take the energy to be about 100 keV, the Vh energy to be about 1-10 MeV, and ~ 
10^^. It is important to emphasize that the growth rate does not depend on the density and 
temperature of a plasma. Then for the growth rate (10) we obtain 7mJl ~ 1.7 — 5.4- lO^sec"^. 
Comparing this with the collision frequency 5 of neutrinos, it can be seen the latter to be 
much less than the former. Indeed, for the given mean free path l^, — 2km{10/e„)'^ — 2-lO^cm 
[7], the coUision frequency is 5 ~ ~ r ~ lO^sec"^, i.e. under these circumstances, the 
medium for the neutrino is " coUisionless" . 
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IV. NONLINEAR WAVES 



We now discuss the nonlinear features of evolution of the unstable waves. If we consider 
the case, when the period of oscillations, is much less than the wave amplitude 

growth time, l/Imfl, then the dispersion relation can be written in a one-dimensional per- 
turbation as — quQi,. This corresponds to the condition for the neglect of the third 
derivatives in Eq.(7). Doing so, we obtain 

d%i 

-^ + u.-u. = 0, (11) 

which describes a beam of noninteracting neutrinos. 

As was shown in [13], the equation type of Eq.(ll) has solutions with the characteristic 
properties of nonlinear waves. It has been verified that the initial perturbation of the velocity 
{u{z, 0) ~ sinkz) in the phase plane {u, z) after some time leads the wave u{z, t) to break. 

Thus, our results indicate that the neutrino beam of noninteracting particles, similar to 
a beam of plasma particles, has many properties of nonlinear system. Namely, it is subject 
to wave breaking and to generation of higher harmonics. 

If we take into account the temperature of neutrinos and suppose that the neutrino 
gas is an ideal {Pn — rinTi, — Poui'ni./no,^)'^), then the equation of motion (7) modifies as 

^ = -^^Vn.-/3V(VV), (12) 

at riv 

where Cl{ni,) — CQ^,{ni,/nou)'^''^, C'oi/ is the velocity of sound in the neutrino gas, and 7 is 
the effective adiabatic index. 

Equations (1) and (12) are so called Boussinesq equations. We can see that the character 
of nonlinear processes strongly depends on the dispersion, i.e., the dependence of phase ve- 
locity on the wave number. The most important nonlinear effect here, as it was shown above, 
is that the wave front steepens. At this point the dispersion becomes significant. Waves 
having different wave number have different phase velocities, and the nonlinear steepening 
of front can be compensated by the dispersive spread. Therefore, in this case the stationary 
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waves can exist, which propagate with constant velocities without changing their shape (e.g., 
sohtons). To note, Eqs.(l) and (12) describe also the nonstationary waves [14]. 

Equations type of Eqs.(l) and (12) have been thoroughly studied, and it has been shown 
that these kind of equations have many periodic as well as solitary solutions. If we use 
the common method for the weak dispersive medium, we obtain from Eqs.(l) and (12) the 
Korteweg-de Vries equation for the neutrino velocity 

du d . ^, 

where u = ^"^^^u^, and a = 2c^- This equation is well studied and we do not repeat all 
solutions and conservation laws of it. Our aim was to prove that the neutrino gas exhibits 
the same features as the waves in plasmas (e.g. electromagnetic waves, waves on the surface 
of a heavy liquid, waves in biology, etc.). 

V. ION DYNAMICS AND INSTABILITIES 

We next consider the excitation of an ion-sound waves by the neutrino flux. Above we 
studied the high-frequency oscillations neglecting the ion dynamics. However, the ion impact 
turns out to be extremely important in the low-frequency waves which may be excited in a 
strong non-isothermal plasma with hot electrons and cold ions. As for the electrons we may 
assume that the electrons are in equilibrium under the conditions of low-frequency oscillation 
and that their density is determined by the Boltzmann formula: 

etj)—y/2Gjr5nu 

He = riQee ^ Srii, ^n^- Uq^ . (14) 

The ions can be described hydrodynamically by means of the HD equations (Tj = 0) : 

= -eV0 , ^ div{niVi) = , (15) 

and the Poisson equation in this case is 

V = 47re(noee - Ui) . (16) 



After linearization of the set of Eqs.(l),(2),(15) and (16) we obtain the dispersion relation 
for the ion-sound oscillation in the presence of the monoenergetic neutrino beam as 



^2^2 , ,2^,^2^2 



where a = is simply the ion-sound frequency, 

QVs /-.ox 
= /, 2 2 ' ^^^) 

here i'. = \/— is the ion-sound velocity. 

We examine the dispersion relation (17) for two interesting cases. Let us first discuss the 
small wavelength case, gru 3> 1 and fl <^ uJs- In this case Eq.(17) reduces to 

{n - quo^f + aq^c^ = , (19) 

which has a solution of the form Q — ReQ, + ilmO, with 

i?efi = quouCosO and ImD, — ^faqc . (20) 

Hence, this solution describes the emission of low-frequency neutrinos inside a resonance 
cone {cosff). 

Next, for the long wavelengths, qrY) <^ 1, the dispersion relation (17) casts in to the form 

(Q2 - q\'':){Sl - qu^r^f - aq\lq\^ = . (21) 

Here we look for the solution with coincide roots, i.e. Q = qVg + F and = qd^^ -|- F, and 
get for the growth rate of the oscillatory instabilities the following expression 

/".r=^(|a)%c. (22) 

VI. FORMATION OF SHOCK WAVES 

Let us now investigate the nonlinear low-frequency waves. To this end, we add the ion 
pressure to the equation of motion of the ions, Eq.(15), and since f2 -C Go's neglect the inertial 
term. Then solution of this equation is 
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Ui = noiC . (23) 



If we further impose the quasi-neutrahty, Ug — rii, we obtain 

ect>^V2GFSn,-^ . (24) 

Substituting (24) into the Boltzmann distribution function of the electrons (14), we get 

— . (25) 

noe 

For the neutrino gas the natural requirement is GpSrii, -C Tg, so that we have 

rie - npe ^ V2GFSn^ ^^g) 
riQe Te + Ti 

This equation when substituted into Eq. (2) yields the equation of motion of the neutrinos 

^ + {u.- V)u. = SV^ , (27) 

where S = ^^Ji^y^T^+r ) ■ ^® linearize this equation and the equation of density (1), then 
we obtain the dispersion equation (19) for Tj = 0. 

We now consider the one-dimensional nonlinear waves, for this case equations are 

— + -„^u. = (28) 
^+-4"' = 4-- (29) 
These equations have exact solutions for the initial profile of the neutrino density 



ch?- 

and in nondimensional quantities are written as 



Mz,0) = ^, (30) 



Q^(l + Q^r^)ch-\C-VT) (31) 
V = -2QTth{^ - Vt) , (32) 

1 /2 

where Q = r = ujpit, ^ = ^z, V = ^, with ujpi = (^^^f^) being the Langmuir 

frequency of ions. 
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Prom Eqs.(31) and (32), we can readily derive 



Examining Eq.(33), we now discuss the question how the maximum of the initial density 
changes in time. For this purpose, we take a derivative of Eq.(33) by Q, i.e. ^ |^=o, and 
allow it to be infinite. Prom this condition, we find the expression for the maximum density 

QmaxiO, ^) = - VT^^) ■ (34) 

One can see that at r ^ 0, Qmax 1, ^1/(0,0) = uqu- Then with increase of time (r) 
the shape of the density changes, namely becomes narrow. Also, the center of the density 
= 0) increases with rise of r and reaches another maximum, Qmax = 2 at t — 1/2. After 
a certain value r — tq, the shock waves are formed, and derivatives are ^ |^=o,t=to= ^ ? 

|^=o,T=To= oo. However, the density itself remains continuous. Using these conditions 
from Eq.(33) follows the relation QqTq = 3/2, where tq > 0.5 and Q^Tq 3> 1. 

3D numerical simulation which underscores the above picture is carried out. Namely, 
we consider 3D problem, when the initial profile of density of neutrinos is the Gaussian 
distribution 

2^1 2^J ■ ^^^^ 

Numerical studies have shown that the dynamics of non-linearity is determined by the 
dimensionality of problem and of the initial profile of the neutrino density. 

In order to confirm the vahdity of ID analytical calculation, we first discuss 3D numerical 
analysis of Eqs.(l) and (27) for pancake-shaped beams, i.e., Oq = 2, and Zq — 0.5. Pigures 
la and lb show that the width of the density becomes shorter in time along the propagation 
direction (z-axis). The density of neutrinos is concentrated in the small region (on the 
z-axis), and at t = 4 starts to form the shock wave. Whereas Pigures Ic and Id show 
the insignificant expansion of the pancake beams in the transverse direction. Therefore, 
we can conclude that the ID analytical calculation is very good approximation for the 3D 
pancake-shaped beams. 
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Quite different processes develop when the initial distribution of neutrinos has a form of 
the bullet, i.e., when qq = 0.5 and z — 2. We can see that there are two special competitive 
regimes: one, which is due to the self-focusing, is the compression of the neutrino density 
in the transverse direction (see Figures 2a and 2b), and the second is the expansion of 
the density along the propagation direction (see Figures 2b and 2c). The redistribution of 
neutrinos in the transverse direction is accompanied by the formation of sharp maximum of 
density near the z-axis (r = \/x^~+y^ — > 0), and as Figure 2a shows, the formation of the 
shock wave takes place at t=2.8. 

VII. SUMMARY AND DISCUSSIONS 

We have studied the problem of weak Fermi interaction of a neutrino beam with a 
plasma. Novel fluid instabilities driven by the neutrino flux in a plasma are observed. 
The range of wavelengths corresponding to several instabilities, discussed in this paper, for 
relevant parameters of the neutrino gas and plasma is A ~ 1 — 10~^cm. We have found that 
for the neutrino beam there exists a Cherenkov type of emission of low energy neutrinos; 
this is different from the usual Cherenkov effect, though. It should be emphasized that in 
our case the physical mechanism is that the high energy neutrino does not just excite the 
plasma wave but also can scatter from density perturbations of electrons and can emit the 
low energy neutrino. We also found that in the usual hot plasma environments such as 
supernova explosions, the effect considered above prevents neutrino escape. Since the low 
energy neutrinos cannot stream away from the central regions of the plasma one can have 
an accumulation of such neutrinos in a dense stellar medium. Thus, the above discussed 
processes can considerably modify certain phases of evolution in a stellar model. It should 
be emphasized that unlike the previous models [11], [16] on neutrino pair emission by a 
stellar plasma, the emission of low energy neutrinos by high energy neutrinos, as discussed 
in this paper, does not depend on the density and temperature of a plasma. The growth 
rate found for some typical neutrino parameters is 5.4 • lO^sec"^. Which is much larger than 
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the collision frequency of neutrinos, lO^sec"^. Thus in this case the medium for the neutrino 
is " coUisionless" . 

The nonlinear stationary as well as nonstationary waves in the neutrino gas are also 
discussed. The dispersive effect, which is responsible for the existence of the nonlinear 
stationary waves, e.g. solitons, in the neutrino gas is due to the electron density modulation 
by the neutrino ponderomotive force [4]. Our results indicate that a neutrino beam is subject 
to wave breaking and shock wave formation. Thus, high energy neutrinos will lose energy also 
by wave breaking in addition to Cherenkov emission, leading to plasma heating (and neutrino 
cooling). Whereas, the shock waves can produce a relativistic energy flow. In addition, the 
solitons can also be a potential candidates for the generation of relativistic particles. These 
particles then release the energy and produce the observed radiation in gamma-ray bursts 
(GRBs). Hence, one has a origin of radiation in the modeling of burst sources, which requires 
a discussion of particle acceleration processes. Since the simplest, most conventional, and 
practically inevitable interpretation of the observations of GRBs is that GRBs result from the 
conversion of the kinetic energy of ultra-relativistic particles to radiation. It is well accepted 
that many of them originate in the very distant, early universe. In the early universe, 
the processes discussed in the paper can also lead to formation of nonlinear structures, 
contributing to the formation of the large scale structure of the Universe [15]. Which, is 
believed, grew gravitationally out of small density fluctuations [17]. The effect of this density 
variation in the early universe was left on the cosmic microwave background radiation in the 
form of spatial temperature fluctuations. It should be emphasized that the gravity, however, 
cannot produce these fluctuations, but increase alone. Therefore a discussion of the physical 
models of generation of the initial matter density fluctuations is very crucial. There is no 
comprehensive model at present that can explain their origin. The weak Fermi interaction in 
a plasma, as discussed in this paper, can be an alternative and a new source for the required 
density perturbations. As demonstrated, a long lived nonlinear structures, which carry large 
amount of mass and energy, arc generated in such a system. Since an initial localization of 
mass and energy is exactly that the gravity needs for eventual structure formation, weak 
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Fermi interaction may have provided a decisive element in the formation of a large scale 
map of the observable Universe. 



APPENDIX A: THE GENERALIZATION OF THE WIGNER-MOYAL 

EQUATION 

The Wigner-Moyal equation in quantum kinetic theory is the analogy of the one-particle 
Liouville equation in kinetic theory In order to derive this equation, we start with a disper- 
sion relation of a single neutrino, taking into account the universal weak Fermi interaction 
between plasma electrons and neutrinos {U — s/^Gprie) 



where E is the energy, p is the momentum, rrii, is the neutrino rest mass. This equation has 
two distinct energy solutions 



The physical meaning of this result is best understood when we are in the rest frame of the 
background medium, which creates the potential energy U. The positive sign in Eq.(A2) 
corresponds to the neutrino solution and it shows that, for a given magnitude of the particle 
momentum, the energy of the neutrino is increased by an amount U with respect to its value 
in vacuum. The other solution in Eq.(A2), corresponding to the negative sign, is that of an 
antineutrino solution. 

As mentioned in the Introduction, in an isotropic plasma one can ignore the spin of 
neutrinos, and therefore we can consider a scalar wave function ■0,^, associated with the 
neutrino gas, and we can write in equilibrium 



{E - Uf = p^c^ + mlc^ , 



(Al) 




(A2) 




(A3) 



where uj 



■|, k = |, /co = and a;^ = ^ is the Fermi frequency. 
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On the other hand, if the neutrino gas is not in equihbrium, the amphtude of ipi, will 
slowly change in space and time due to the interaction of neutrinos with the background 
plasma fluctuations. We can then use the geometrical optic approximation, a; — > a; — id/dt, 
k ^ k + iV and Ue — noe + (uj ^\ \, \ k |^| V |). Doing expansion 



/ ^ N 

V + ^ V)2 + ki^Y.i )i^^ + ^lf~' (2^^V - V')^' (A4) 



N 
j=0 

and neglecting higher order derivatives greater than in Eq.(A4), we arrive from Eq.(A3) 
to 

<| + «.V)^„ - |(V^ - (|V)>„ - a,,^VV = , (A5) 

where use was made of the relation uj — c{k'^ + k^Y^"^ + Upinoe) , and the following notations 
were introduced 

' {k^ + A;^)V2 ' I dn ^n=noe ~ noe ' ^ ^ 

We now write the product of two wave functions which are defined at two distinct points 

and instants of time 

F = Mri, ti)ruir2, t2) = ^.(^ + ^, t + ^)^: - ^, t - ^) , (A7) 

where we have introduced new space and time variables defined by 

-R = ^(ri + r2), f=ri-f2, t^^{ti+t2), t ^ h - t2 . 

Let us make a Fourier transform of (A7) on the variables r, r in order to introduce the 
Wigner function f{R,t,k,uj) 

f{R,t,k,uj) = j^J drj dTi;.{R+ ^)^:(^- ^je-^^^^-— ) . (A8) 

This is the generalized Wigner distribution function. Integration of (A8) by the frequency 
UJ leads to the ordinary expression of the Wigner distribution function [18] as 



e-'"' . (A9) 
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Prom (A8) also follows important relation 

Assuming that the neutrinos verify the dispersion relation ( A2) , we can say that for a given 
value k we always have a well-defined value of the frequency, uj — uj{k). This means that we 
can write 

f{R, t, k, uj) = 27r f{R, t, k)S{uj - uj{k)) . (All) 
Integration over the wave vector spectrum leads to the density of neutrinos 

— * 

n.{R,t) =1 ij. \'= J ^J{R,t,k) . (A12) 

Following the procedure described in previous works [18-21,4], we can derive an equation 
for the function (A7), the result is 

i{| + v,Vn}F - |{V« • V. - (|V^) (|v.)}f - ^(5n, - 6n,)F = . (A13) 

Now making the Fourier transformation of Eq.(A13), we obtain an evolution equation for 
function (AS) in the form 

{- + vl"WR)f{R,t,k,uj) j dfj driSn, - 5n2)Fe-*('^™) = , (A14) 

where 



ifff = ^ n J -0 



1+ ^' 
" ^ "0 



' k^ + kl 



If all derivatives of 5n exist, we can expand Sni{R + |,t + |) — Sn2{R — — ^) around 
{R,t),i.e. 

5m -Sm = pj{'^-Vn + r|)'[l - (-l)']5n(^, t) (A15) 



and obtain the generalized Wigner-Moyal equation 
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Integration of Eq.(Al6) by the frequency leads to the Wigner-Moyal equation for the Wigner 
function (A9) without the last term on the right hand side of Eq.(Al6). 

If we just retain the first term on the right hand side of Eq.(Al6), i.e., the term corre- 
sponding to sinx ~ X, we obtain the well-known Liouville-Vlasov kinetic equation for the 
distribution function (A8), but with an additional time derivative term 

(I . 5/'V«)/(«. , - ;^{v«.n . V./ - ^ . f ) = . (A17) 

The stationary solution of Eq.(Al7) is the Fermi distribution function 



fp = (e + l) . (A18) 

Prom Eq.(Al7) we can now derive a set of fluid equations for the neutrino gas. To this 
end we shall introduce the principal definitions of quantities which describe the state of the 
neutrinos in a plasma. Namely, the total number of neutrinos is defined as 



N, 



— * 

= J J 7^3 J ^fiR,t,k,u;) = J dRn.{R,t) , (A19) 



(27r) 

the total energy of the neutrino gas 



— * 

Etotai - JdRj I ^nujf{R,t,k,u;) = I dRs{R,t) , (A20) 

and the neutrino mean velocity 

I f dk r duj kc^ ^, -> . , , . 

J [2t{Y J ztt uj 

Having definitions (A19-A21) and using the usual momentum procedure, from Eq.(Al7) we 
obtain a set of fluid equations 



dt 



+ div{n^u^) = , (A22) 



— + K-V)k^ = -- -(V + — , A23 

at < ujk > ^ ot^ noe n„ 



^ + V • = hivpn,-^ , A24 
at at rioe 
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where we used the definitions of the momentum density 



= ^ / ^/ ^kf{R,t,k,u;) , (A25) 



and 

rij, f dk f duj f{R,t,k,u!) 



<ujk> J (27r)3 y 27r a; ■ ^ ^ 

The quantity P^, plays role of the neutrino pressure and is defined as 

Equations (A22-A24) together with equations of the plasma, given in this paper, describe 
the collective process, which plays a more important role than the pair collision effect. 
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FIGURES 

FIG. 1. a) and b) The time evolution of the neutrino density along the z - axis in the case of 
the pancake-shaped beams 

c) and d) The time evolution of the neutrino density along the x - axis in the case of the pan- 
cake-shaped beams 

FIG. 2. a) The time evolution of the neutrino density along the x - axis in the case of the 
neutrino bullet 

b) The neutrino density at t=2.8 in (z,x) space 

c) The time evolution of the neutrino density along the z - axis in the case of the neutrino bullet 
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Fig. la Tsintsadze 




Fig. lb Tsintsadze 




Fig.lc Tsintsadze 




Fig. Id Tsintsadze 




Fig. 2a Tsintsadze 




Fig. 2b Tsintsadze 




Fig. 2c Tsintsadze 



